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We have performed a search for bursts of gravitational waves associated with the very bright Gamma Ray 
Burst GRB030329, using the two detectors at the LIGO Hanford Observatory. Our search covered the most sen- 
sitive frequency range of the LIGO detectors (approximately 80-2048 Hz), and we specifically targeted signals 
shorter than ~150 ms. Our search algorithm looks for excess correlated power between the two interferome- 
ters and thus makes minimal assumptions about the gravitational waveform. We observed no candidates with 
gravitational wave signal strength larger than a pre-determined threshold. We report frequency dependent upper 
limits on the strength of the gravitational waves associated with GRB030329. Near the most sensitive frequency 
region, around ~250 Hz, our root-sum-square (RSS) gravitational wave strain sensitivity for optimally polarized 
bursts was better than h_Rss— 6x 10^^^ Hz"^''^. Our result is comparable to the best published results searching 
for association between gravitational waves and GRBs. 

PACS numbers: 04.80.Nn, 07.05.Kf, 95.85.Sz 97.60.Bw 
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I. INTRODUCTION 

Gamma Ray Bursts (GRBs) are short but very energetic 
pulses of gamma rays from astrophysical sources, with du- 
ration ranging between 10 ms and 100 s. GRBs are histori- 
cally divided into two classes (11 0] based on their duration: 
"short" (< 2 s) and "long" (> 2 s). Both classes are isotropi- 
cally distributed and their detection rate can be as large as one 
event per day. The present consensus is that long GRBs Q| 
are the result of the core collapse of massive stars resulting 
in black hole formation. The violent formation of black holes 
has long been proposed as a potential source of gravitational 
waves. Therefore, we have reason to expect strong association 
between GRBs and gravitational waves 0,|3,|5|. In this pa- 
per, we report on a search for a possible short burst of gravita- 
tional waves associated with GRB030329 using data collected 
by the Laser Interferometer Gravitational Wave Observatory 
(LIGO). 

On March 29, 2003, instrumentation aboard the HETE- 
2 satellite 0] detected a very bright GRB, designated 
GRB030329. The GRB was followed by a bright and well- 
measured afterglow from which a redshift Q] of z =0.1685 
(distance ~800 Mpc |0]) was determined. After approxi- 
mately 10 days, the afterglow faded to reveal an underlying 
supernova (SN) spectrum, SN2003dh This GRB is the 
best studied to date, and confirms the link between long GRBs 
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and supernovae. 

At the time of GRB030329, LIGO was engaged in a 2- 
month long data run. The LIGO detector array consists of 
three interferometers, two at the Hanford, WA site and one at 
the Livingston, LA site. Unfortunately, the Livingston inter- 
ferometer was not operating at the time of the GRB; therefore, 
the results presented here are based on the data from only the 
two Hanford interferometers. The LIGO detectors are still 
undergoing commissioning, but at the time of GRB030329, 
their sensitivity over the frequency band 80 to 2048 Hz ex- 
ceeded that of any previous gravitational wave search, with 
the lowest strain noise 

of -6x10^22 Hz-i/2 around 250 Hz. 
A number of long GRBs have been associated with X-ray, 
radio and/or optical afterglows, and the cosmological origin 
of the host galaxies of their afterglows has been unambigu- 
ously established by their observed redshifts, which are of or- 
der unity 01 ■ The smallest observed redshift of an optical af- 
terglow associated with a detected GRB (GRB980425 flolfTll 
[ij) is z=0.0085 (^35 Mpc). GRB emissions are very likely 
strongly beamed ll3Lfl4tl . a factor that affects estimates of the 
energy released in gamma rays (a few times 10^° erg), and 
their local true event rate (about 1 per year within a distance 
of lOOMpc). 

In this search, we have chosen to look for a burst of 
gravitational waves in a model independent way. Core col- 
lapse 01, black hole formation fisll and black hole ring- 
down Il6l [ItIi may each produce gravitational wave emis- 
sions, but there are no accurate or comprehensive predictions 
describing the gravitational wave signals that might be asso- 
ciated with GRB t ype sources. Thus, a traditional matched 
filtering approach Il8l [l9il is not possible in this case. To 
circumvent the uncertainties in the waveforms, our algorithm 
does not presume any detailed knowledge of the gravitational 
waveform and we only apply general bounds on the waveform 
parameters. Based on current theoretical considerations, we 
anticipate the signals in our detectors to be weak, comparable 
to or less than the detector's noise j20ll2lll22il . 

This paper is organized as follows: Section II summarizes 
the currently favored theories of GRBs and their consequences 
for gravitational wave detection. Section III provides obser- 
vational details pertinent to GRB030329. Section IV briefly 
describes the LIGO detectors and their data. Section V dis- 
cusses the method of analysis of the LIGO data. In Section VI 
we compare the events in the signal region with expectations 
and we use simulated signal waveforms to determine detec- 
tion efficiencies. We also present and interpret the results in 
this section. Section VII offers a comparison with previous 
analyses, a conclusion, and an outlook for future searches of 
this type. 



II. PRODUCTION OF GRAVITATIONAL WAVES IN 
MASSIVE CORE COLLAPSES 

The apparent spatial association of GRB afterglows with 
spiral arms, and by implication star formation regions in re- 
mote galaxies, has lead to the current "collapsar" or "hyper- 
nova" scenario l23l l24ll in which the collapse of a rotating. 



4 



massive star to a Kerr black hole can lead to relativistic ejecta 
emitted along a rotation axis and the associated production 
of a GRB jet. The identification of GRB030329 with the su- 
pernova SN2003dh (section 3 below) gives further support to 
this association. This observation is consistent with the the- 
ory that the GRB itself is produced by an ultra-relativistic jet 
associated with a central black hole. Stellar mass black holes 
in supernovae must come from more massive stars. Ref. ll24ll 
presents "maps" in the metallicity-progenitor mass plane of 
the end-states of stellar evolution and shows that progenitors 
with 25 Mq can produce black holes by fall-back accretion. 

The observed pulsar kick velocities of ~500 km/s hint at 
a strong asymmetry around the time of maximum compres- 
sion, which may indicate deviations from spherical symme- 
try in the progenitor. The resulting back reaction on the core 
from the neutrino heating provides yet another potential phys- 
ical mechanism for generating a gravitational wave signal. In 
the model of |25] it imparts a kick of 400-600 km/s and an 
induced gravitational wave strain roughly an order of magni- 
tude l arg er than in j 20i| and an order of magnitude smaller 
than II2I. 

Theoretical work on gravitational wave (GW) signals in the 
process of core-collapse in massive stars has advanced much 
in recent years, but still does not provide detailed waveforms. 
Current models take advantage of the increase in computa- 
tional power and more sophisticated input physics to include 
both 2D and 3D calculations, utilizing realistic pre-collapse 
core models and a detailed, complex equation of state of su- 
pernovae that produce neutron stars. The most recent studies 
by independent groups give predictions for the strain ampli- 
tude within a similar range, despite the fact that the domi- 
nant physical mechanisms for gravitational wave emission in 
these studies are different t20. 21. 26. 27, 281. The calcu- 
lations of I20I1 are qualitatively different from previous core 
collapse simulations in that the dominant contribution to the 
gravitational wave signal is neutrino-driven convection, about 
20 times larger than the axisymmetric core bounce gravita- 
tional wave signal. 

The applicability of the above models to GRBs is not 
clear, since the model endpoints are generally neutron stars, 
rather than black holes. Another recent model involves ac- 
cretion disks around Kerr black holes |2§|, subject to non- 
axisymmetric Papaloizou-Pringle instabilities I30ll in which 
an acoustic wave propagates toroidally within the fallback 
material. They are very interesting since they predict much 
higher amplitudes for the gravitational wave emission. 

For our search, the main conclusion to draw is that in spite 
of the dramatic improvement in the theoretical models, there 
are no gravitational waveforms that could be reliably used as 
templates for a matched filter search, and that any search for 
gravitational waves should ideally be as waveform indepen- 
dent as practical. Conversely, detection of gravitational waves 
associated with a GRB would almost certainly provide cru- 
cial new input for GRB/SN astrophysics. It is also clear that 
the predictions of gravitational wave amplitudes are uncertain 
by several orders of magnitude, making it difficult to predict 
the probability to observe the gravitational wave signature of 
distant GRBs. 
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Time from GRB030329 Fregate_B trigger (T„ = 732973047.7) [s] 

FIG. 1: The GRB030329 light curve as measured by the HETE-2 
FREGATE B instrument. The arrow indicates the HETE trigger time. 
The signal region analyzed in this study is indicated by the horizontal 
bar at the top. This figure is the courtesy of the HETE collaboration. 

The timeliness of searching for a gravitational wave signal 
associated with GRBs is keen in light of the recent work by 
1I21I1 and l2(ill . (20'] finds that the signal due to neutrino con- 
vection exceeds that due to the core bounce and therefore a 
chaotic signal would be expected. Studies with simplified or 
no neutrino transport (e.g., I21I1 . ll22il ') find the core-bounce 
to be the dominant contributor to the GW signal. The large- 
scale, coherent mass motions involved in the core bounce 
leads to a predicted gravitational wave signal resembling a 
damped sinusoid. 



III. GRB030329 RELATED OBSERVATIONAL RESULTS 
A. Discovery of GRB 030329 and its afterglow 

On March 29, 2003 at 11:37:14.67 UTC, a GRB trig- 
gered the FREGATE instrument on board the HETE-2 satel- 
lite 13 m m The GRB had an effective duration of 
~50 s, and a fluence of 1 .08 x lO^"^ erg/cm^ in the 30-400 keV 
band JH. The KONUS detector on board the Wind satel- 
lite also detected it t34i. triggering about 15 seconds after 
HETE-2. KONUS observed the GRB for about 35 seconds, 
and measured a fluence of 1.6x10^'* erg / cm^ in the 15- 
5000 keV band. The measured gamma ray fluences place this 
burst among the brightest GRBs. Figure^shows the HETE-2 
light curve for GRB030329 IH- 

The rapid localization of the GRB by HETE ground analy- 
sis gave an accurate position which was distributed about 73 
minutes after the original trigger. A few hours later, an opti- 
cal afterglow [2, 36] was discovered with magnitude R=12.4, 
making it the brightest optical counterpart to any GRB de- 
tected to date. The RXTE iIstIi satellite measured a X-ray 
flux of 1.4x 10"^" erg s"^ cm~^ in the 2-10 keV band about 
4h51m after the HETE trigger, making this one of the bright- 
est X-ray afterglows detected by RXTE 1:38]. The National 
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Radio Astronomy Observatory (NRAO) observed the 
radio afterglow, which was the brightest radio afterglow de- 
tected to date 14011. Sp ectroscopic measurements of the bright 
optical afterglow j4lll revealed emission and absorption lines, 
and the inferred redshift [z = 0.1685, luminosity distance 
Dl ^ 800 Mpc) made this the second nearest GRB with 
a measured distance. To date, no host galaxy has been identi- 
fied. It is likely that numerous other GRBs have been closer 
than GRB030329, but the lack of identified optical counter- 
parts has left their distances undetermined. 

Spectroscopic measurements I^ l4^l43|l . about a week after 
the GRB trigger, revealed evidence of a supernova spectrum 
emerging from the light of the bright optical afterglow, which 
was designated SN2003dh. The emerging supernova spec- 
trum was similar to the spectrum of SN1998bw a week before 
its brightness maximum L44tii45t1 . 

SN1998bw was a supernova that has been spatially and 
temporally associated with GRB980425 E El III, and was 
located in a spiral arm of the barred spiral galaxy ESO 184- 
G82 at a redshift of z = 0.0085 (~35 Mpc), making it the 
nearest GRB with a measured distance. The observed spec- 
tra of SN2003dh and SN1998bw, with their lack of hydro- 
gen and helium features, place them in the Type Ic super- 
nova class. These observations, together with the observa- 
tions linking GRB980425 (which had a duration of ~23 s) 
to SN1998bw, make the case that collapsars are progenitors 
for long GRBs more convincing. In the case of SN1998bw, 
Woosley et al. j46ll and Iwamoto et al. ifTlll found that its ob- 
served optical properties can be well modeled by the core col- 
lapse of a Ch-O core of mass 6 M0 (main sequence mass of 
25 Mq) with a kinetic energy of ~2x 10^^ ergs. This energy 
release is about an order of magnitude larger than the energies 
associated with typical supernovae. 



to a value larger than Qj. Effectively, prior to this time the 
relativistic ejecta appears to be part of a spherical expansion, 
the edges of which cannot be seen because they are outside of 
the beam, while after this time the observer perceives a jet of 
finite width. 

This leads to a faster decline in the light curves. Zeh et al. 
and Li et al. j48l l49ll show that the initial "break" or strong 
steepness in the light curve occurs at about 10 hours after the 
initial HETE-2 detection. 

Frail et al. lfl3ll give a parametric relation between beaming 
angle <dj, break time tj, and Ei^o as: 

^ \24hours J \ 2 J 

^ yiO^^ergs) V0.2/ V0.1cto-3/ 

where Qj is measured in radians. It was argued that the fire- 
ball converts the energy in the ejecta into gamma rays effi- 
ciently llsoll (r]j saO.2), and that the mean circumburst den- 
sity is nwO.l cm~3 15111 . Evaluating equation 13.21 for the 
parameters of GRB030329 (tj w 10 hours, z=0.1685, and 
Eiso=2x 10^2 erg) gives 9^ w0.07 rad. 

Therefore the beaming factor that relates the actual energy 
released in gamma rays (E-^) to the isotropic equivalent en- 
ergy is 9|/2si 1/400, so that E^ w5xlO'*^ erg. Compa ring 
Eiso and with the histograms in Fig. 2 of Frail et al. 113 
, GRB030329 resides at the lower end of the energy distribu- 
tions. The calculated isotropic energy from GRB980425, the 
GRB associated with SN1998bw, is also low (~10''® erg). 



IV. OVERVIEW OF THE LIGO DETECTORS 



B. GRB030329 energetics 

A widely used albeit naive quantity to describe the energy 
emitted by GRBs is the total isotropic equivalent energy in 
gamma rays: 

E^so - MBC)Dlf/il + z) w 2 X 10^2 (3 J) 

where / is the measured fluence in the HETE-2 waveband 
and BC is the approximate bolometric correction for HETE-2 
for long GRBs. Using a "Band spectrum" iITtIi with a single 
power law to model the gamma ray spectrum, and using a 
spectral index, /3 = —2.5, gives that the GRB's total energy 
integrated from 1 keV to 5 GeV is greater than that present in 
the band 30-400 keV by a factor 2.2. 

However, it is generally believed that GRBs are strongly 
beamed, and that the change in slope in the afterglow light 
curve corresponds to the time when enough deceleration has 
occurred so that relativistic beaming is diminished to the point 
at which we "see" the edge of the jet. This occurs during 
the time in which the relativistic ejecta associated with the 
GRB plows through the interstellar medium, and the beaming 
factor r^^, where F is the bulk Lorentz factor of the flow, 
increases from a value smaller than the beaming angle Qj, 



The three LIGO detectors are orthogonal arm Michelson 
laser interferometers, aiming to detect gravitational waves by 
interferometrically monitoring the relative (differential) sep- 
aration of mirrors, which play the role of test masses. The 
LIGO Hanford Observatory (LHO) operates two identically 
oriented interferometric detectors, which share a common 
vacuum envelope: one having 4 km long arms (HI), and one 
having 2 km long arms (H2). The LIGO Livingston Observa- 
tory operates a single 4 km long detector (LI). The two sites 
are separated by ~3000 km, representing a maximum arrival 
time difference of ~±10 ms. 

A complete description of the LIGO interferometers as they 
were configured during LIGO's first Science Run (SI) can be 
found in Ref IH. 



A. Detector calibration and configuration 

To calibrate the error signal, the response to a known dif- 
ferential arm strain is measured, and the frequency-dependent 
effect of the feedback loop gain is measured and compensated 
for During detector operation, changes in calibration are 
tracked by injecting continuous, fixed-amplitude sinusoidal 
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FIG. 2: Typical LIGO Hanford sensitivity curves during the S2 Run 
[strain Hz^^''^] (black and grey lines). The LIGO design sensitivity 
goals (SRD) are also indicated (dotted and dashed lines). 



excitations into the end test mass control systems, and moni- 
toring the ampUtude of these signals at the measurement (er- 
ror) point. Calibration uncertainties at the Hanford detectors 
were estimated to be < 11%. 

Significant improvements were made to the LIGO detectors 
following the S 1 run, held in early fall of 2002: 

1. The analog suspension controllers on the H2 and LI 
interferometers were replaced with digital suspension 
controllers of the type installed on HI during SI, re- 
sulting in lower electronics noise. 

2. The noise from the optical lever servo that damps the 
angular excitations of the interferometer optics was re- 
duced. 

3. The wavefront sensing system for the HI interferometer 
was used to control 8 of 10 alignment degrees of free- 
dom for the main interferometer As a result, it main- 
tained a much more uniform operating point over the 
run. 

4. The high frequency sensitivity was improved by oper- 
ating the interferometers with higher effective power, 
about 1.5 W. 

These changes led to a significant improvement in detec- 
tor sensitivity. Figure |2] shows typical spectra achieved by 
the LIGO interferometers during the S2 run. The differences 
among the three LIGO spectra reflect differences in the oper- 
ating parameters and hardware implementations of the three 
instruments which are in various stages of reaching the final 
design configuration. 



B. The second science run 

The data analyzed in this paper were taken during LIGO's 
second Science Run (S2), which spanned approximately 60 
days from February 14 to April 14, 2003. During this time, 
operators and scientific monitors attempted to maintain con- 
tinuous low noise operation. The duty cycle for the interfer- 
ometers, defined as the fraction of the total run time when 
the interferometer was locked and in its low noise configura- 
tion, was approximately 74% for HI and 58% for H2. The 
longest continuous locked stretch for any interferometer dur- 
ing S2 was 66 hours for HI . 

At the time of the GRB030329 both Hanford interferom- 
eters were locked and taking science mode data. For this 
analysis we relied on the single, ~4.5 hours long coincident 
lock stretch, which started ~3.5 hours before the trigger time. 
With the exception of the signal region, we utilized ~98% of 
the data within this lock stretch as the background region (de- 
fined in section V). 60 seconds of data before and after the sig- 
nal region were not included in the background region. Data 
from the beginning and from the end of the lock stretch were 
not included in the background region to avoid using possi- 
bly non-stationary data, which might be associated with these 
regions. 

As described below, the false alarm rate estimate, based on 
background data, must be applicable to the data within the 
signal region. We made a conservative choice and avoided 
using background data outside of the lock stretch containing 
the GRB trigger time. This is important when considering the 
present non-stationary behavior of the interferometric detec- 
tors. 



V. ANALYSIS 

The goal of the analysis is either to identify significant 
events in the signal region or, in the absence of significant 
events, to set a limit on the strength of the associated gravita- 
tional wave signal. Simulations and background data were 
used to determine the detection efficiency for various ad- 
hoc and model-based waveforms (Section fVIBt and the false 
alarm rate of the detection algorithm respectively. 

The analysis takes advantage of the information provided 
by the astrophysical trigger The trigger time determined 
when to perform the analysis. As discussed below, the time 
window to be analyzed around the trigger time was chosen to 
accommodate most current theoretical predictions and timing 
uncertainties. The source direction was needed to calculate 
the attenuation due to the LIGO detector's antenna pattern for 
the astrophysical interpretation. 

The two co-located and co-aligned Hanford detectors had 
very similar frequency-dependent response functions at the 
time of the trigger Consequently, the detected aiTival time and 
recorded waveforms of a gravitational wave signal should be 
essentially the same in both detectors. It is natural then to con- 
sider cross-correlation of the two data streams as the basis of 
a search algorithm. This conclusion can also be reached via a 
more formal argument based on the maximum log-Likelihood 
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The schematic of the full analysis pipeline is shown in Fig- 
ure |3l The underlying analysis algorithm is described in de- 
tail in Ref. l54ll . The background data, the signal region data 
and the simulations are all processed identically. The back- 
ground region consists of the data where we do not expect 
to have a gravitational wave signal associated with the GRB. 
We scan the background to determine the false alarm distribu- 
tion and to set a threshold on the event strength that will yield 
an acceptable false alarm rate. This threshold is used when 
scanning the signal region and simulations. In order to esti- 
mate our sensitivity to gravitational waves, simulated signals 
of varying strength are added to the detector data streams. 

The signal region around the GRB trigger is scanned to 
identify outstanding signals. If events were detected above 
threshold, in this region, their properties would be tested 
against those expected from gravitational waves. If no events 
were found above threshold, we would use the estimated sen- 
sitivity to set an upper limit on the gravitational wave strain at 
the detector 

The output from each interferometer is divided into 330 sec 
long segments with a 15 second overlap between consecu- 
tive segments (both ends), providing a tiling of the data with 
300 second long segments. In order to avoid edge effects, 
the 180 sec long signal region lies in the middle of one such 
300 sec long segment. This tiling method also allows for adap- 
tive data conditioning and places the conditioning filter (see 
Sec. IV B II below) transients well outside of the 300 second 
long segment containing the signal region. 



A. Choice of signal region 

Current models suggest 0] that the gravitational wave sig- 
nature should appear close to the GRB trigger time. We con- 
servatively chose the duration and position of the signal region 
to over-cover most predictions and to allow for the expected 
uncertainties associated with the GRB trigger timing. A 180 
second long window (see Figure 0, starting 120 seconds be- 



fore the GRB trigger time is sufficient; roughly ten times 
wider than the GRB light curve features, and wide enough to 
include most astrophysical predictions. Most models favor an 
ordering where the arrival of the gravitational wave precedes 
the GRB trigger (3], but in a few other cases the gravitational 
wave arrival is predicted to be contemporaneous Ia. l5^ to the 
arrival and duration of the gamma rays (i.e after the GRB trig- 
ger). The 60 second region after the GRB trigger time, is suf- 
ficient to cover these predictions and also contains allowance 
for up to 30 seconds uncertainty on trigger timing, which is a 
reasonable choice in the context of the HETE light curve. Fig- 
ure shows a signal rise time of order 10 s, precursor signals 
separated from the main peak, and significant structure within 
the main signal itself. Effects due to the beaming dynamics 
of the GRB and the instrumental definition of the trigger time 
can also be significant contributors to the timing uncertainty. 



B. Search algorithm 

1. Data Conditioning 

The data-conditioning step was designed to remove instru- 
mental artifacts from the data streams. We used an identical 
data conditioning procedure when processing the background, 
the signal region and the simulations. 

The raw data streams have narrowband lines, associated 
with the power line harmonics at multiples of 60Hz, the vi- 
olin modes of the mirror suspension wires and other narrow 
band noise sources. The presence of lines has a detrimental 
effect on our sensitivity because lines can produce spurious 
correlations between detectors. In addition, the broadband 
noise shows significant variations over timescales of hours 
and smaller variations over timescales of minutes and seconds 
due to alignment drift and fluctuations. The background data 
must portray a representative sample of the detector behavior 
around the time of the trigger. Broadband non-stationarity can 
limit the duration of this useful background data and hence the 
reliability of our estimated false alarm rate. 

Our cross correlation based algorithm performs best on 
white spectra without line features. We use notch filters to 
remove the well-known lines, such as power line and vio- 
lin mode harmonics from both data streams. Strong lines of 
unknown origin with stationary mean frequency are also re- 
moved at this point. We also apply a small correction to miti- 
gate the difference between the phase and amplitude response 
of the two Hanford detectors. 

We bandpass filter and decimate the data to a sampling rate 
of 4096 Hz to restrict the frequency content to the ~80 Hz to 
~2048 Hz region, which was the most sensitive band for both 
LIGO Hanford detectors during the S2 run. 

In order to properly remove weaker stationary lines and the 
small residuals of notched strong lines, correct for small slow 
changes in the spectral sensitivity and whiten the spectrum of 
the data, we use adaptive line removal and whitening. As all 
strong lines are removed before the adaptive whitening, we 
avoid potential problems due to non-stationary lines and en- 
hance the efficiency of the follow up adaptive filtering stage. 
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The conditioned data has a consistent white spectrum without 
major lines and sufficient stationarity, from segment to seg- 
ment, throughout the background and signal regions. 

The end result of the pre-processing is a data segment with 
a flat power spectral density (white noise), between ~80 Hz 
and ~2048 Hz. The data conditioning was applied consis- 
tently after the signal injections. This ensures that any change 
in detection efficiency due to the pre-processing is properly 
taken into account. 



2. Gravitational Wave Search Algorithm 

The test statistics for a pair of data streams are constructed 
as follows. We take pairs of short segments, one from each 
stretch, and compute their cross-correlation function. The ac- 
tual form of the cross-correlation used (Cr''" ) is identical to 
the common Euclidean inner product: 

_ j 

^Kv,3 ^ ^™ [fc + [A: + p + , (5.1) 



where the pre-conditioned time series from detector "x" is 
Hj; = {iJ^[0], . . .} and i,k,p and j are all integers in- 

dexing the data time series, with each datum being (1/4096) s 
long. As we now only consider the two Hanford detectors "m" 
and "n" can only assume values of 1 (i?i) or 2 {H^). There 
are therefore three free parameters to scan when searching for 
coherent segments of data between a pair of interferometers 
(m,n); 1 . the center time of the segment from the first detector 
(k); 2. the relative time lag between the segments from the 
two detectors (p); and 3. the common duration of segments 
(2jH-l) called the integration length. 

The optimum integration length to use for computing the 
cross-correlation depends on the duration of the signal and 
its signal-to-noise ratio, neither of which is a priori known. 
Therefore the cross-correlation should be computed from seg- 
ment pairs with start times and lengths varying over values, 
which should, respectively, cover the expected arrival times 
(signal region) and consider durations of the gravitational 
wave burst signals 1 20. ,2 1 ., ,25, ,26. .27., 2M( --0(l-128ms)). 

Hence we apply a search algorithm i54ll that processes the 
data in the following way. 

(1) A three dimensional quantity {Cu.j [p]) is constructed: 



1/2 



(5.2) 



scanning the range of segment center times (k), integration 
lengths (2jH-l) and relative time shifts (p = 0,±l,±2,...).A 
coherent and coincident signal is expected to leave its local- 
ized signature within this three dimensional quantity. 

We use a fine rectangular grid in relative time shift (p) and 
integration length (2jH-l) space. The spacing between grid 
points is ~1 ms for the segment center time (k) and (1/4096) s 
for the relative time shift. The spacing of the integration 
lengths is approximately logarithmic. Each consecutive inte- 
gration length is ~50% longer than the previous one, covering 
integration lengths from ~1 ms to ~128 ms. 



Introducing small, non-physical relative time shifts (much 
larger than the expected signal duration) between the two data 
streams before computing the cross-correlation matrix sup- 
presses the average contribution from a GW signal. This prop- 
erty can be used to estimate the local noise properties, thereby 
mitigating the effects of non- stationarity in the interferometer 
outputs. Accordingly, Cfe j \p\ contains the autocorrelation of 
the coherent signal for relative time shifts at and near p = 
(called "core"), while far away, in the "side lobes", the contri- 
bution from the signal autocorrelation is absent, sampling only 
the random contributions to the cross-correlation arising from 
the noise. The optimal choice of the core size depends on the 
expected signal duration (integration length), the underlying 
detector noise and it cannot be smaller than the relative phase 
uncertainty of the datastreams. The core region can reach as 
far as 5 ms, as it increases with increasing integration length. 
The size of each side lobe is twice the size of the core region 
and the median time shift associated with the side lobes can 
be as large as 325 ms as it is also increasing with increasing 
integration length. We use the side lobes of Ck,j [p] to esti- 
mate the mean {p,k.j) and variance {dk.j) of the local noise 
distribution, which is also useful in countering the effects of 
non-stationarity. 

(2) The three dimensional quantity is reduced to a two di- 
mensional image (see Fig.|4}, called a corrgmm, as follows. 
The values of Ck.j [p] in the core region are standardized by 
subtracting Jikj and then dividing by 'Skj- Positive standard- 
ized values in the core region are summed over p to determine 
the value of the corrgram pixel. Each pixel is a measure of the 
excess cross-correlation in the core region when compared to 
the expected distribution characterized by the side lobes for 
the given (k,j) combination. 

(3) A list of events is found by recursively identifying and 
characterizing significant regions (called "clusters") in the 
corrgram image. Each event is described by its arrival time, its 
optimal integration length and its strength (ES). The event's 
arrival time and its optimal integration length correspond to 
the most significant pixel of the cluster. The event strength 
is determined by averaging the five most significant pixels of 
the cluster, as this is helpful in discriminating against random 
fluctuations of the background noise. 

The strength of each event is then compared to a preset de- 
tection threshold corresponding to the desired false alarm rate. 
This detection threshold is determined via extensive scans of 
the background region. 



VI. RESULTS 
A. False alarm rate measurements 

In order to assess the significance of the cross-correlated 
power of an event, we determined the false alarm rate versus 
event strength distribution. We used the full background data 
stretch for this measurement. 

Figure |5] shows the event rate as a function of the event 
strength threshold for the background region. The error bars 
reflect 90% CL Poisson errors, based on the the number of 
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FIG. 4: Examples of corrgram images. The horizontal axes are time 
(linearly scaled) and the vertical axes are integration length (loga- 
rithmically scaled). The color axis, an indicator of the excess corre- 
lation, is independently auto-scaled for each quadrant for better vis- 
ibility, therefore the meaning of colors differ from quadrant to quad- 
rant. The time ticks also change from quadrant to quadrant for better 
visibility. The rainbow type color scale goes from blue to red, dark 
red marking the most significant points within a quadrant. The upper 
two quadrants show the corrgram image of injected Sine-Gaussians 
(250 Hz, Q = 8.9). The bottom quadrants are examples of noise. The 
maximum of the intensity scale is significantly higher for both quad- 
rants with injections, when compared to the noise examples. The top 
left injection is strong enough to be significantly above the preset de- 
tection threshold, while the top right injection is weak enough to fall 
significantly below the detection threshold. 



events within the given bin. We used this distribution to fix 
the event strength threshold used in the subsequent analysis. 

We chose an event strength threshold with an associated 
false alarm rate of less than ~5x 10~^ Hz, equivalent to less 
than ~9% chance for a false alarm within the 180 second long 
signal region. 



B. Efficiency determination 

The detection sensitivity of the analysis was determined 
by simultaneously adding simulated signals of various ampli- 
tudes and waveforms to both data streams in the background 
region and evaluating the efficiency of their detection as a 
function of the injected amplitude and waveform type. 

The waveforms we considered include Sine-Gaussians to 
emulate short narrow -band bursts, Gaussians to emulate short 
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FIG. 5: False alarm rate as a function of the event strength threshold 
as determined from background data. The error bars reflect 90% CL 
Poisson errors, based on the the number of events within the given 
bin. The pointer indicates the event strength threshold used for the 
analysis, which corresponds to an interpolated false alarm rate of less 
than 5 x 10^'* Hz. Note that the signal region data is not included in 
this calculation. The position of symbols correspond to the center of 
the bins. 

broad-band signals, and Dimmelmeier-Font-Miiller numerical 
waveforms 12611 . as examples of astrophysically motivated sig- 
nals. 

Calibration of the waveforms from strain to ADC counts 
was performed in the frequency domain, and was done sepa- 
rately for each interferometer. Calibration procedures of the 
LIGO detectors are described in Ref. j52il . The transformed 
signals, now in units of counts of raw interferometer noise, 
were then simply added to the raw data stream. 

The amplitudes and the times of the injections were ran- 
domly varied. In this way we ensured that each amplitude 
region sampled the full, representative range of noise varia- 
tions and that we had no systematic effects, for example, due 
to a regular spacing in time. 

To a reasonable approximation the sensitivity of our anal- 
ysis pipeline can be expressed in terms of the frequency con- 
tent, the duration and the strength of the gravitational wave 
signal. Therefore, it is sufficient to estimate the sensitivity of 
our search for a representative set of broad and narrow band 
waveforms, which span the range of frequencies, bandwidth, 
and duration we wish to search. 

We characterize the strength of an arbitrary waveform by its 
root-sum-square amplitude (/i_r55), which is defined as 15611 : 



hRss = ^J \h{t)\^dt. (6.1) 

The above definition of hfj^ss includes all frequencies, while 
the gravitational wave detectors and search algorithms are 
only sensitive in a restricted frequency band. In principle, 
one can analogously define a "band-limited" h^ss, in which 
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FIG. 6: Efficiency of the detection algorithm for a sample wave- 
form as a function of signal strength (huss)', in this case a Sine- 
Gaussian of fo = 250 Hz and Q = 8.9. To extract this curve numer- 
ous simulated waveforms were embedded in a representative fraction 
of the background data at random times with randomly varying sig- 
nal strength. The plot shows the fraction of signals detected as the 
function of amplitude and a sigmoid function fit. The reconstructed 
signal onset times were required to fall within ±60 ms of the true 
onsets, which also explains why the low hiiss end of the curve falls 
near zero. This is a typical plot and in general, the agreement be- 
tween the measured values, and the fit is better than ~5%. We relied 
on the fit to extract our upper limits for an optimally oriented and 
polarized source. Section lVI El below describes the corrections due 
to non-optimal source direction and polarization. 

only the sensitive frequency band of the analysis is taken into 
account. Within this paper we choose to adopt the Eq. 16.11 
definition of h^jss for historical reasons. 

The extracted sensitivities (see the examples in Figures |6] 
and can be used to generalize our measurements and es- 
timate the pipeline's sensitivity for other similar band-Utnited 
waveforms. 

To assess the sensitivity for relatively narrow band wave- 
forms, we used Sine-Gaussian injections of the form: 

h{t) = ho sin(cjo i)e-*'/2'"' , (6.2) 

with a central angular frequency of uiq = 27r/o, and Q = 
u!o(J = 27r/ocr. The relation between ho and hj^ss is given 
as; 



The injected signals covered the frequency range between 
100 and 1850 Hz with 13 values of fo- To test the depen- 
dence of the sensitivity on signal duration, we used three val- 
ues of Q (4.5, 8.9 and 18) for each frequency (see Table 



FIG. 7: Detected event strength versus Hrss of the injected Sine- 
Gaussian waveform with fo = 250 Hz and Q ~ 8.9. The dots in- 
dicate the scatter of the distribution of raw measurements. The gray 
band shows the quadratic polynomial fit, which allows us to convert 
the strength of an observed event into the equivalent hnss value and 
determine the associated 90% CL error bars. The markers with er- 
ror bars represent the 90% CL regions for subsets of the data. For 
each marker, 90% of the measurements used were within the hori- 
zontal error bars and 90% of the detected event strengths values fell 
within the vertical error bars. The vertical dash-dot line represents 
the 50% detection efficiency associated with the waveform type and 
the chosen detection threshold (horizontal dotted line). As expected, 
the crossing of the threshold and the 50% efficiency lines agree well 
with the fit and the center of the corresponding marker. The verti- 
cal dashed line represents the boundary of the region where we have 
better than 90% detection efficiency. The "comer" defined by the 
event strength threshold and the 90% detection efficiency boundary 
(dashed lines) agrees well with the curve outlined by the lower end 
of the vertical error bars of the markers. All events in the upper right 
corner of the plot (above and beyond the dashed lines) are detectable 
with high confidence. This plot is typical for different waveforms 
considered in the analysis. 



Near the most sensitive frequency region, around ~250 Hz, 
our gravitational wave strain sensitivity for optimally polar- 
ized bursts was better than h/{55~5 X 10~^^ Hz~^/^. Figure|8l 
shows the sensitivity for these narrow band waveforms. The 
symbols mark the simulated event strength (hjiss) necessary 
to achieve 90% detection efficiency for each waveform. We 
quote the gravitational wave signal strength associated with 
the 90% detection efficiency, as this can be related to the up- 
per limits on the gravitational wave strength associated with 
the source. Figure |S] also illustrates the insensitivity of the 
detection efficiency to the Q of the Sine-Gaussian waveforms 
with the same central frequency, as these reach their 90% ef- 
ficiency levels at similar gravitational wave strengths, even 
though their Q differ by a factor of ~4; for a given hj^ss, 
a longer signal (higher Q) would of course, have a smaller 
hpEAK- This strength is frequency dependent, naturally fol- 
lowing the frequency dependence of the detector sensitivities, 
which are also indicated in Figure|8] 

Table im shows a similar set of efficiencies estimated using 
broad-band simulated signals. We used two types of broad- 
band waveforms, Sine-Gaussians with unity quality factor and 
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FIG. 8: Sensitivity of the detection algorithm for detecting Sine- 
Gaussian waveforms versus characteristic frequency. The plot shows 
the strength necessary for 90% detection efficiency. The grey spectra 
illustrate the sensitivity of the 2K and 4K Hanford detectors during 
the time surrounding the GRB030329 trigger. The error bars reflect 
a total 15% error. 



Gaussians. Both are short bursts, however, the Gaussians are 
even functions while the Sine-Gaussians are odd, leading to 
different peak amplitudes with the same huss value. Gaus- 
sians were parametrized as: 

h{t) = /loe-*'/'"' . (6.4) 
The relationship between ho and hfjss for a Gaussian is: 
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FIG. 9: Number of events versus event strength in the signal re- 
gion (circle). The diamonds show the expected distribution based 
on the background region. The squares mark the expected distribu- 
tion based on non-physical time shifts (ranging from 2 to 9 seconds) 
between the HI and H2 datastreams in the background region. The 
error bars reflect 90% CL Poisson errors. The position of the sym- 
bols correspond to the center of the bins. 



placed an upper limit on the detected strength of gravitational 
waves associated to GRB030329. Our fixed false alarm rate 
permitted the results of simulations to be used directly in set- 
ting upper limits. 

The upper limits on h fiss for narrow band waveforms are 
given in TableU TablesHllandlllllshow the upper limits for the 
broadband simulations and astrophysically motivated wave- 
forms, respectively. 



D. Errors 



h%ss = ho y/v^ . (6.5) 

The estimated sensitivities indicate that the 90% detection 
efficiency limits for short bursts are similar to those obtained 
for the narrow band waveforms when one takes into account 
that only part of the power of the broad-band waveforms is 
confined to the analysis frequency band. Longer Gaussian 
bursts are more difficult to detect, as their spectrum has a sig- 
nificant low frequency component, outside the sensitive band 
of our analysis. 

We have also estimated our efficiency for a set of astrophys- 
ically motivated burst waveforms lE6ll (see Table IIIH . These 
simulated waveforms are not expected to be necessarily asso- 
ciated with GRBs, rather these results are presented here to 
further illustrate the waveform independence of the analysis. 

C. Signal region 

The analysis of the signal region (Fig. |9} yielded only 
events well below the predetermined event strength threshold 
(<60% of threshold). Since we had no candidate event, we 



The analysis method, the procedures used to determine the 
efficiencies, and the non-stationary nature of the data, all con- 
tribute to the uncertainty associated with the results. 

The efficiency (versus /i/fss) values have an estimated 
~11% uncertainty due to our limited knowledge of the cal- 
ibrated response of our detectors. This estimate also accounts 
for the slight difference in calibrated response between the 
signal region and background data used for the simulations. 

An additional uncertainty arises from the non-stationarity 
of the data. The results of the simulations exhibit a slight de- 
pendence on the choice of the actual data segments ("base" 
data) used for the injections. This dependence was character- 
ized via simulations using numerous different sub-segments 
of the background data. We repeated the full efficiency es- 
timation process several times for the same waveform, while 
injecting into various base data stretches. The variation in the 
measured upper limits indicated ~10% uncertainty due to the 
dependency of our upper limits on the base data. This un- 
certainty shall also account for the statistical error due to the 
finite number of simulations used. 

We characterized the detection efficiencies for each wave- 
form considered via fits of sigmoid functions (see for example 
Figure |6|l. The fits agree well with the data, but small differ- 
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ences are occasionally observed in the > 90% efficiency re- 
gion. We estimate that using these fits can underestimate the 
90% hmits by < 5%. 

The uncertainties listed above are taken into account by 
specifying a total 15% uncertainty for each measurement in 
Figure|8]and in all Tables. 

The false alarm rate associated with the results was also 
measured. The false alarm rate limit is based on the mea- 
surement with zero lag data plus the 90% confidence Poisson 
error bars. We have checked the assumption of Poisson back- 
ground statistics by examining the time intervals between con- 
secutive triggers and the variance in trigger counts for varying 
ES thresholds when the background sample is divided into 
50 equal-length intervals. Good agreement with the Poisson 
expectation is observed. This choice provides a conservative 
estimate of our associated (~5 x 10"'* Hz) false alarm rate. 

E. Astrophysical interpretation 

GRB030329 has a well-determined redshift, therefore we 
can relate our observed limits on strain to a measure of the 
total gravitational wave energy emission. For a strain h{t) at 
distance Dj^ from a source of gravitational radiation, the as- 
sociated power is proportional to (h = dh/dt), though the 
proportionality constant will depend on the (unknown) emis- 
sion pattern of the source and the antenna pattern of the detec- 
tor (for the known source position, but unknown polarization 
angle). 

In general, it is not possible to relate our upper limit on 
the strain from a particular waveform to a Umit on the energy 
radiated by the source, without assuming a model. Sources 
that radiate energy i?GW might produce an arbitrarily small 
signal h{t) in the detector, e.g., if the dynamics in the source 
were purely axi-symmetric with the detector located on the 
axis. Nevertheless, we can associate a strain h{t) in the de- 
tector with some minimum amount of gravitational-wave en- 
ergy radiated by the source by choosing an "optimistic" emis- 
sion pattern, thereby obtaining a measure of the minimum 
amount of energy that would need to be radiated in order to 
obtain a detectable signal. We will show that the progenitor 
of GRB030329 is not expected to have produced a detectable 
signal. 

We are interested in a "plausible case scenario" of gravita- 
tional wave emission in order to obtain the minimum (plausi- 
ble) amount of gravitational-wave energy radiated that could 
be associated with a detector signal h{t). We do not expect the 
gravitational waves to be strongly beamed, and we expect that 
we are observing the GRB progenitor along some preferred 
axis. We take a model best case scenario to be that of gravi- 
tational wave emission from a triaxial ellipsoid rotating about 
the same axis as the GRB (i.e., the direction to the Earth). If 
we assume quadrupolar gravitational wave emission, the plus- 
and cross-polarization waveforms, emitted at a polar angle 9 
from the axis of rotation to be: 

h+ = i(l + 008^61) /i+,o (6.6) 
hx = cosOhx.o (6.7) 



where /i+.o and h^.o are two orthogonal waveforms (e.g., a 
Sine-Gaussian and a Cosine-Gaussian), each containing the 
same amount of radiative power That is, we assume that the 
same amount of gravitational-wave energy is carried in the 
two polarizations and that they are orthogonal: 

h^ Qdt— / h\ f^dt and / h+fihx,odt ^Q. 

-OO J — OO J —OO 

(6.8) 

Thus, we would expect that the gravitational waves travelling 
along the rotational axis (toward the Earth) would be circu- 
larly polarized, and that the detector would receive the signal 

h^F+h+,o + Fxhx,n (6.9) 

where and represent the detector responses to the po- 
larization components ft,+.o ™d hx,o IstIi . and depend on the 
position of the source in the sky and on a polarization angle. 
The radiated energy from such a system is calculated to be 

(6.10) 

where = + F^ (which depends only on the position 
of the source on the sky) and where we are integrating over a 
spherical shell around the source with radius Dl (the distance 
to the Earth). Alternatively, using Parseval's identity, we have 

o 2 3 r)2 foc 

Egw^^^^ \fh\'df (6.11) 

where 

/OO 
h{t)e-^'"^* dt. (6.12) 
-OO 

Whereas optimal orientation gives ?; = 1 for a source at 
zenith, the position of GRB030329 was far from optimal. The 
angle with respect to zenith was 68° and the azimuth with re- 
spect to the a;-arm was 45°, which yields rj ~ 0.37. 

We now relate E'gw to the strain upper limits using the 
specific waveforms used in the analysis. For a Gaussian wave- 
form [see Eq. ( I6.4H : 

and for a sine-Gaussian waveform [see Eq. ( 16. 2H : 

---(#)(^)<'-«^-"'' 

where Q = oJo(J = 271/0(7. The relation between ho and 
hfiss is given in Eqs. ( 16. 5t and ( 16. 3t . 

We can relate the observed limit on Hbss to an equivalent 
mass A/eq which is converted to gravitational radiation with 
100% efficiency, Sqw = J^Ieqc^, at a luminosity distance 
Dl ~ 800 Mpc. For sine-Gaussian waveforms with fo = 
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250 Hz and Q = 8.9, Meq = 1.9 x W^tj-^Mq. For Gaus- 
sian waveforms with a = 1ms, Meq = 3.1 x W^ti^^AIq. 
However, we would not expect that the gravitational-wave lu- 
minosity of the source could exceed ~ c^/G = 2 x IO^A/qc^ 
per second [58], so we would not expect an energy in gravita- 
tional waves much more than ~ 2 x 10'^ Mqc^ in the ~ 10 ms 
Sine-Gaussian waveform, or an energy of much more than 
~ 3 X IQ'^Mqc^ in the maximum duration (150 ms) of the 
search; far below the limits on Meqc^ that we find in this 
analysis. Present theoretical expectations on the gravitational 
wave energy emitted range from 10-6 Mq c2 - 10"^ Mq c^ to 
10~^ M0 c^ - MfT^ c^ for some of the most optimistic models 
[see e.g. ISL Ba. l59l 16011 1. Nonetheless, these scalings indi- 
cate how we can probe well below these energetic limits with 
future analyses. For example, assuming similar detector per- 
formance for an optimally oriented trigger like GRB980425 
(Dl w35 Mpc) the limit on the equivalent mass would be 
^^EQ ~60 M0 for the Gaussian waveforms mentioned above 
with cr = 1 ms. 



VII. SUMMARY 

A. Comparison with previous searches for gravitational waves 
from GRBs 

Our result is comparable to the best published results 
searching for association between gravitational waves and 
GRBs 16111 . however these studies differ in their most sensitive 
frequency. 

Tricarico et al. l63l used a single resonant mass detector, 
AURIGA l63ll , to look for an excess in coincidences be- 
tween the arrival times of GRBs in the BATSE 4B catalog. 
They used two different methods. They searched for events 
identified above a certain threshold in the gravitational wave 
data, and also attempted to establish a statistical association 
between GRBs and gravitational waves. No significant excess 
was found with the former method. The latter used a vari- 
ant of the correlation based Finn-Mohanty-Romano (FMR) 
method l64ll . However, instead of using the cross-correlation 
of two detectors, as proposed in the FMR method, only the 
variance of the single detector output was used. A sample of 
variances from times when there were no GRBs was com- 
pared with a corresponding sample from data that spanned 
the arrival times of the GRBs. An upper limit on the source- 
averaged gravitational wave signal root mean square value of 
1.5 X lO-^'* was found using 120 GRBs. This limit applies at 
the AURIGA resonant frequencies of 913 and 931 Hz, which 
are very far from the most sensitive frequency of the LIGO de- 
tectors (~250 Hz). This work l62il was later extended 16^ . 
which led to an improved upper limit. 

The data analysis method employed in Modestino & Mo- 
leti 16^ is another variant of the FMR method. Instead of 
constructing off-source samples from data segments that are 
far removed from the GRB trigger, the off-source samples 
are constructed by introducing non-zero time shifts between 
the two detector data streams and computing their cross- 
correlation. For narrowband resonant mass detectors, the di- 



rectional information of a GRB cannot be exploited to dis- 
criminate against incorrect relative timing since the signal in 
the output of the detector is spread out by the detector re- 
sponse over time scales much larger than the light travel time 
between the detectors. 

Astone et al. l67l l6^ report on a search for a statistical 
association between GRBs and gravitational waves using data 
from the resonant mass detectors EXPLORER j 69] and NAU- 
TILUS |7(|]. They report a Bayesian upper limit on gravita- 
tional wave signal amplitudes of 1.2x10^^^, at 95% proba- 
bility, when the maximum delay between the GRB and grav- 
itational wave is kept at 400 sec. The upper limit improves 
to 6.5 X 10"^^ when the maximum delay is reduced to 5 sec. 
However, the absence of directional and/or distance informa- 
tion for most of these GRBs precluded accounting for source 
variations; the gravitational wave signal amplitude was as- 
sumed to be the same for all of the GRBs. 

Astone et al. jFjlll report on the operation of the resonant 
mass detectors EXPLORER during the closest ever gamma 
ray burst (GRB980425) with known redshift and direction. At 
the time of the burst, EXPLORER was taking data with close 
to optimal orientation. GRB980425 was ~23 times closer to 
Earth than GRB030329 giving a ~520 increase in energy sen- 
sitivity. Based on their sensitivity and the loudest event within 
±5 minutes of the GRB980425 trigger the authors quote a 
limit of ~1600 M© for a simple model assuming isotropic 
gravitational wave emission. 

Recently, Astone et al. |6 1 1 executed a search aiming to 
detect a statistical association between the GRBs detected by 
the satellite experiments BATSE and BeppoSAX, and the EX- 
PLORER and NAUTILUS gravitational wave detectors. No 
association was uncovered. Their upper limit is the lowest 
published result, which is based on bar-detector gravitational 
wave data. 



B. Conclusion 

We have executed a cross-correlation-based search for pos- 
sible gravitational wave signatures around the GRB030329 
trigger, which occurred during the Second Science Run (S2) 
of the LIGO detectors. We analyzed a 180 second signal re- 
gion around the GRB and 4.5 hours of background data, sur- 
rounding the signal region, corresponding to a single coinci- 
dent lock stretch. These data were sufficient to characterize 
the background, scan the signal region and estimate our effi- 
ciency. We used the same procedure, based on cross correla- 
tion, for each of these studies. We evaluated the sensitivity of 
the search to a large number of broad and narrow band wave- 
forms. 

We observed no candidates with gravitational wave signal 
strength larger than a pre-determined threshold, therefore we 
set upper limits on the associated gravitational wave strength 
at the detectors. The present analysis covers the most sensitive 
frequency range of the Hanford detectors, approximately from 
80 Hz to 2048 Hz. The frequency dependent sensitivity of our 
search was hi^55~0(6x lO'^i) Hz-^/^. 

The prospect for future searches is promising, as the sen- 
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sitivity of the instruments improves with further commission- 
ing. 

Once operating at target sensitivity, the detectors will be 
more sensitive to strain than they were during S2 by factors of 
10 - 100, depending on frequency (see Figure|2]). This implies 
an improvement of a factor of ~1000 in sensitivity to Egw, 
since Eqw scales like ^ hj^ss (^^^ for example Eg. 16.131 1. 

Detection of GRBs with measured redshifts signifi- 
cantly smaller than GRB030329's is certainly possible. 
GRB030329's electromagnetic brightness was due to a favor- 
able combination of distance and our position in its beam. One 
year of observation will incorporate hundreds of GRBs with 
LIGO data cov erag e and some of these GRBs, even though 
fainter IH IHItI than GRB030329, could be significantly 
closer, as was 1998bw. We can also hope for sources with 
more optimal direction and coincidence between three or four 
observing interferometers. 

VIII. ACKNOWLEDGEMENTS 

The authors gratefully acknowledge the support of the 
United States National Science Foundation for the construc- 



tion and operation of the LIGO Laboratory and the Particle 
Physics and Astronomy Research Council of the United King- 
dom, the Max-Planck-Society and the State of Niedersach- 
sen/Germany for support of the construction and operation of 
the GEO600 detector The authors also gratefully acknowl- 
edge the support of the research by these agencies and by the 
Australian Research Council, the Natural Sciences and Engi- 
neering Research Council of Canada, the Council of Scientific 
and Industrial Research of India, the Department of Science 
and Technology of India, the Spanish Ministerio de Ciencia 
y Tecnologia, the John Simon Guggenheim Foundation, the 
Leverhulme Trust, the David and Lucile Packard Foundation, 
the Research Corporation, and the Alfred P. Sloan Founda- 
tion. We are grateful to Scott Barthelmy and the GCN net- 
work and Kevin Hurley and the IPN network for providing us 
with near real time GRB triggers and to the Ulysses, Konus, 
SAX, and HETE experiments, who detect and generate the 
events distributed by GCN and IPN. This research has made 
use of data obtained from the HETE science team via the web- 
site http://space.mit.edu/HETE/Bursts/Data. HETE is an in- 
ternational mission of the NASA Explorer program, run by 
the Massachusetts Institute of Technology. 



[1] C. Kouveliotou, et al., Astrophysical Journal, Letters 413, LlOl 
(1993). 

[2] B. Zhang and P. Meszai'os, Int. J. Mod. Phys. A (2003), in press, 
astro-ph/03 11321. 

[3] C. L. Fryer, S. E. Woosley, and A. Heger, Astrophysical Journal 
550, 372 (2001). 

[4] M. B. Davies, et al., Astrophysical Journal 579, L63 (2002). 

[5] M. H. van Putten, et al.. Physical Review D 69 (2004). 

[6] D. Q. Lamb, et al.. New Astronomy Review 48, 423 (2004). 

[7] P. A. Price, et al.. Nature 423, 844 (2003). 

[8] T. Matheson, et al., Astrophysical Journal 599, 394 (2003). 

[9] V. V. Sokolov, et al., ArXiv Astrophysics e-prints (2003), astro- 
ph/03 12359. 

[10] S. R. Kulkami, et al.. Nature 395, 663 (1998). 
[1 1] K. Iwamoto, et al.. Nature 395, 672 (1998). 
[12] T. J. Galama, et al.. Nature 395, 670 (1998). 
[13] D. A. Frail, et al., Astrophysical Journal, Letters 562, L55 
(2001). 

[14] M. H. P. M. van Putten and T. Regimbau, Astrophysical Journal, 

Letters 593, L15 (2003). 
[15] K. Belczynski, T. Bulik, and B. Rudak, Astrophysical Journal, 

Letters 608, L45 (2004). 
[16] E. E. Flanagan and S. A. Hughes, Physical Review D 57, 4535 

(1998). 

[17] J. D. E. Creighton, Physical Review D 60 (1999). 

[18] C. W. Helstrom, Statistical Theory of Signal Detection (Perga- 

mon Press, London, 1968), 2nd ed. 
[19] B. Abbott, et al.. Physical Review D 69 (2004). 
[20] E. Muller, et al.. The Astrophysical Journal 603, 221 (2004). 
[21] C. L. Fryer and M. S. Warren, The Astrophysical Journal 601, 

391 (2004). 

[22] C. D. Ott, et al.. The Astrophysical Journal 600, 834 (2004). 
[23] A. I. MacFadyen and S. E. Woosley, The Astrophysical Journal 
524, 262 (1999). 



[24] A. Heger, et al., The Astrophysical Journal 591, 288 (2003). 
[25] A. Burrows and J. Hayes, Physical Review Letters 76, 352 

(1996) . 

[26] H. Dimmelmeier, J. A. Font, and E. Miiller, Astronomy & As- 
trophysics 393, 523 (2002). 

[27] C. L. Fryer, D. E. Holz, and S. A. Hughes, The Astrophysical 
Journal 565, 430 (2002). 

[28] T. Zwerger and E. Miiller, Astronomy & Astrophysics 320 

(1997) . 

[29] M. H. P M. van Putten, Astrophysical Journal 583, 374 (2003). 
[30] J. C. B. Papaloizou and J. E. Pringle, Monthly notices of the 

Royal Astronomical Society 208, 721 (1984). 
[31] R. Vanderspek, et al., GRB Circular Network 1997 (2003). 
[32] G. R. Ricker, R. K. Vanderspek, and HETE Science Team, 

American Astronomical Society Meeting 203 (2003). 
[33] R. Vanderspek, et al., ArXiv Astrophysics e-prints (2004), 

astro-ph/0401311. 
[34] S. Golenetskii, et al., GRB Circular Network 2026 (2003). 
[35] Http://space.mit.edu/HETE/Bursts/GRB030329. 
[36] Y. M. Lipkin, et al., Astrophysical Journal 606, 381 (2004). 
[37] H. V. Bradt, R. E. Rothschild, and J. H. Swank, Astronomy and 

Astrophysics, Supplement 97, 355 (1993). 
[38] F Marshall and J. H. Swank, GRB Circular Network 1996 

(2003). 

[39] B. E., A. M. Soderberg, and D. A. Frail, GCN notice 2014 
(2003), also see http://www.nrao.edu/pr/2003/grb030329/. 

[40] E. Berger, et al., Astrophysical Journal 599, 408 (2003). 

[41] J. Greiner, et al.. Nature 426, 157 (2003). 

[42] K. Z. Stanek, et al., American Astronomical Society Meeting 
203 (2003). 

[43] P. A. Mazzali, et al., Astrophysical Journal, Letters 599, L95 
(2003). 

[44] J. Hjorth, et al.. Nature 423, 847 (2003). 
[45] J. Greiner, et al. Nature 426, 157 (2003). 



15 



[46] S. E. Woosley, R. G. Eastman, and B. P. Schmidt, Astrophysical 

Journal 516, 788 (1999). 
[47] D. Band, et al., Astrophysical Journal 413, 281 (1993). 
[48] A. Zeh, et al., GRB Circular Networlc 2047 (2003). 
[49] W. Li, et al., GRB Circular Networlc 2064 (2003). 
[50] A. M. Beloborodov, Astrophysical Journal, Letters 539, L25 

(2000) . 

[51] D. A. Frail, E. Waxman, and S. R. Kulkarni, Astrophysical Jour- 
nal 537, 191 (2000). 

[52] LIGO Scientific Collaboration: B. Abbott et al., Nucl. Instrum. 
Methods A 517 (2004). 

[53] W. G. Anderson, et al.. Physical Review D 63 (2001). 

[54] S. Mohanty et al., 8th GWDAW conference proceedings. Class. 
Quantum Grav. (2004). 

[55] R. A. Araya-Gochez, Monthly Notices of the RAS 355, 336 
(2004). 

[56] B. Abbott, et al.. Physical Review D 69 (2004). 

[57] S.W. Hawking and W. Israel, editors, 300 Years of Gravitation 

(Cambridge University Press, 1987), Chapters by K.S. Thome. 
[58] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation 

(Freeman, 1973). 
[59] H. J. Mosquera Cuesta, Physical Review D 65, 064009 (2002). 
[60] R. Ruffini, et al., Astrophysical Journal, Letters 555, LI 13 

(2001) . 

[61] P. Astone, et al., ArXiv Astrophysics e-prints (2004), astro- 



ph/0408544. 

[62] P Tricarico, et al.. Physical Review D 63 (2001). 

[63] G. A. Prodi et al.. Initial operation of the gravitational wave 
detector AURIGA (1998), in the proceedings of the "Sec- 
ond Edoardo Amaldi conference on gravitational wave experi- 
ments", (CERN - Switzerland, 1997), edited by E. Coccia, G. 
Veneziano, G. Pizzella, World Scientific, Singapore p. 148-158. 

[64] L. S. Finn, S. D. Mohanty, and J. Romano, Physical Review D 
60(1999). 

[65] P. Tricarico, A. Ortolan, and P. Fortini, Classical and Quantum 

Gravity 20, 3523 (2003). 
[66] G. Modestino and A. Moleti, Physical Review D 65 (2002). 
[67] P. Astone, et al.. Astronomy and Astrophysics, Supplement 138, 

603 (1999). 

[68] P Astone, et al.. Physical Review D 66 (2002). 
[69] P Astone, et al.. Physical Review D 47, 362 (1993). 
[70] P. Astone, et al., Astroparticle Physics 7, 231 (1997). 
[71] L. Amati, et al.. Astronomy and Astrophysics, Supplement 138, 
605 (1999). 

[72] S. Woosley, Nature 430, 623 (2004), news and views. 

[73] S. Y. Sazonov, A. A. Lutovinov, and R. A. Sunyaev, Nature 
430, 646 (2004), letters to nature. 

[74] A. M. Soderberg, et al.. Nature 430, 648 (2004), letters to na- 
ture. 



16 



TABLE I: hass [Hz~^^'^] for 90% detection efficiency for Sine- 
Gaussians (SG) waveforms at various frequencies (/o) and Q (see 
eg. I6.2> . The quoted values are the results of simulations and are 
subject to a total of ~ 15% statistical and systematic errors, which 
are taken into account when quoting the U Lf^^'J^'" values. Note that 
at the low and at the high frequency end, the low Q waveforms have 
significant power outside of the analysis frequency band. 



Waveform 


a [ms] 


Q 


fa [HZ\ 


1,90% r rr, -l/2l 


UL^^^^ [Hz \ 


SG 


7.2 


4.5 
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17 X lO"^'^ 
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SG 
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4.5 
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SG 


2 
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SG 
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O A 1A — 21 
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554 
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n 1 ^ , in — 21 

9.1 X 10 


SG 
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TABLE II: As in Table|l| hB.ss [Hz'^^^] for 90% detection effi- 
ciency for Gaussian (GA) waveforms of various durations (a) (see 
eQ. I6.4t and for Sine-Gaussians (SG) waveforms at various frequen- 
cies (/o) and Q — 1 (see ea. l6.2L Note that these broadband wave- 
forms have significant power outside of the analysis frequency band. 



Waveform 


a [ms] Q /o [Hz] 
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9.7 X 10"^^ 


SG 


0.35 1 


458 


10 X 10" 


-21 


12 X 10"^^ 


SG 


0.29 1 


554 


13 X 10" 


-21 


14 X 10"^^ 


SG 


0.23 1 


702 


18 X 10" 


-21 


20 X 10-21 


SG 


0.19 1 


850 


23 X 10" 


-21 


26 X 10-21 


SG 


0.16 1 


1000 


26 X 10" 


-21 


30 X 10-21 


SG 


0.12 1 


1361 


39 X 10" 


-21 


45 X 10-21 


SG 


0.11 1 


1458 


44 X 10" 


-21 


51 X 10-21 


SG 


0.1 1 


1554 


46 X lo- 


-21 


52 X 10-21 


SG 


0.094 1 


1702 


ss X 10" 


-21 


63 X 10-21 


SG 


0.086 1 


1850 


61 X 10" 


-21 


70 X 10-21 


GA 


0.5 




8.3 X 10- 


-21 


9.6 X 10-21 


GA 


0.75 




9.6 X 10" 


-21 


1.1 X 10-2° 


GA 


1 




1.3 X 10" 


-20 


1.5 X 10-2° 


GA 


2 




3.3 X 10" 


-20 


3.8 X 10-2° 


GA 


3 




8.2 X 10" 


-20 


9.5 X 10-2° 


GA 


4 




1.9 X 10" 


-19 


2.2 X 10-" 


GA 


5.5 




8.S X 10" 


-19 


9.8 X IQ-i^ 


GA 


8 




1.3 X 10" 


-17 


1.5 X IQ-i^ 


GA 


10 




1.0 X 10" 


-16 


1.2 X IQ-i^ 
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TABLE III: As in Table|l| Hrss [Hz'^^'^] for 90% detection effi- 
ciency for astrophysically motivated waveforms. These waveforms 
are described in detail in Ref. l26ll . Note that most of these wave- 
forms have significant power outside of the analysis frequency band. 

Simulation Waveform h%fs[Hz-^^^] ULl^^l^^jHz-^'^] 



DFM 


AlBlGl 


12 


X 


10" 


-21 


14 


X 


10" 


-21 


DFM 


A1B2G1 


13 


X 


10" 


-21 


15 


X 


10" 


-21 


DFM 


A1B3G1 


12 


X 


10" 


-21 


14 


X 


10" 


21 


DFM 


A1B3G2 


12 


X 


10" 


-21 


14 


X 


10" 


-21 


DFM 


A1B3G3 


12 


X 


10" 


-21 


14 


X 


10" 


-21 


DFM 


A1B3G5 


34 


X 


10" 


-21 


39 


X 


10" 


21 


DFM 


A2B4G1 


24 


X 


10" 


-21 


27 


X 


10" 


-21 


DFM 


A3B1G1 


19 


X 


10" 


-21 


21 


X 


10" 


-21 


DFM 


A3B2G1 


20 


X 


10" 


-21 


23 


X 


10" 


-21 


DFM 


A3B2G2 


15 


X 


10" 


-21 


17 


X 


10" 


-21 


DFM 


A3B2G4 


14 


X 


10" 


-21 


16 


X 


10" 


-21 


DFM 


A3B3G1 


28 


X 


10" 


-21 


33 


X 


10" 


-21 


DFM 


A3B3G2 


17 


X 


10" 


-21 


20 


X 


10" 


-21 


DFM 


A3B3G3 


12 


X 


10" 


-21 


14 


X 


10" 


-21 


DFM 


A3B3G5 


30 


X 


10" 


-21 


34 


X 


10" 


-21 


DFM 


A3B4G2 


23 


X 


10" 


-21 


27 


X 


10" 


-21 


DFM 


A3B5G4 


26 


X 


10" 


-21 


29 


X 


10" 


-21 


DFM 


A4B1G1 


38 


X 


10" 


-21 


44 


X 


10" 


-21 


DFM 


A4B1G2 


32 


X 


10" 


-21 


36 


X 


10" 


-21 


DFM 


A4B2G2 


42 


X 


10" 


-21 


48 


X 


10" 


-21 


DFM 


A4B2G3 


39 


X 


10" 


-21 


45 


X 


10" 


-21 


DFM 


A4B4G4 


17 


X 


10" 


-21 


19 


X 


10" 


-21 


DFM 


A4B4G5 


12 


X 


10" 


-21 


13 


X 


10" 


-21 


DFM 


A4B5G4 


21 


X 


10" 


-21 


25 


X 


10" 


-21 


DFM 


A4B5G5 


19 


X 


10" 


-21 


22 


X 


10" 


-21 



